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ABSTRACT 

A  method  for  utilizing  point  position  radar  data  to  determine  the 
aerodynamic  drag  of  a  projectile  is  described.  Proof  of  the  method's 
validity  and  feasibility  is  represented  by  results  obtained  with  flight 
test  data  taken  for  a  range  firing  of  the  155mm  howitzer  with  the  M107 
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I.  INTRODUCTION 


The  differential  equations  fcr  the  motion  of  a  projectile  are  the  pri¬ 
mary  tools  used  in  the  preparation  of  Firing  Tables.  These  equations  re¬ 
quire  two  basic  kinds  of  input  data:  aerodynamic  coefficients  (determined 
from  aerodynamic  spark  ranges  or  wind  tunnels)  and  ballis;'  t  data  (obtained 
from  large  scale  range  firings).  The  possibility  of  obtaining  both  the  aero¬ 
dynamic  coefficients  and  the  ballistic  data  from  a  single  large  scale  range 
firing  has  beeo  discussed  by  those  who  prepare  Firing  Tables  ever  since 
the  availability  of  tracking  radars  on  many  of  our  national  firing  ranges. 

Of  importance  in  utilising  radar  data  to  obtain  information  for  the  pre¬ 
paration  of  Firing  Tables  are  the  determination  of  suitable  radars  and  the 
development  of  suitable  analysis  methods.  The  Ballistic  Research 
Laboratories  has  been  investigating  these  areas  using  data  supplied  by 
the  following  radars  at  Wallops  Island,  Virginia:  Spandar,  FPS  16, 

MPS  19  and  MOO  XL  For  detailed  descriptions  of  these  radars  see 

a 

Reference  1, 

This  report  describes  a  method  of  utilising  point  position  radar  data 
in  the  determination  of  aerodynamic  drag  ,  the  most  important  of  the 
aerodynamic  coefficients  used  in  the  preparation  of  Firing  Tables.  Also 
included  are  the  results  obtained  when  this  method  was  applied  to  FPS  16 
radar  data  from  a  Bring  of  the  1 55mm  M10?  projectile. 

IL  TEST  PROGRAM 

A  test  program  totaling  28  rounds,  distributed  as  shown,  was  fired 
at  Wallops  Island  to  provide  radar  tracking  data  for  a  series  of  trajec¬ 
tories  which  included  a  typical  sample  of  launch  angles  and  initial 

•References  appear  on  page  31. 

Prectin;  we  Unk 


velocities.  Two  projectiles  were  fired  at  each  of  nine  charge  and  quadrant 
elevation  combinations.  One  projectile  at  each  combination  was  fired  on 
one  day  while  the  second  projectile  at  each  combination  was  fired  on 
another  day  to  obtain  weapon,  atmospheric  and  radar  day  to  day  variations. 
In  addition,  a  single  group  of  10  rounds  was  fired,  as  rapidly  as  possible, 
at  charge  7  and  700  mils  quadrant  elevation  to  provide  a  large  amount  of 
data  on  a  single  occasion  under  similar  meteorological  and  test  conditions. 

Distribution  of  Rounds  Fired 


Quadrant  Elevation-mils 

Charge  Muzzle  Velocity 
(No. )  (mps) 

2  235 

5  372 

7  564 


300  700  1120 

2  2  2 

2  2  2 

2  2,10  2 


Four  radars,  the  FPS  16,  MPS  19,  Spandar  and  MOD  II,  were  used  to 
track  the  test  rounds,  producing  time-space  histories  of  each  projectile’s 
flight.  Additional  data  recorded  for  each  round  were  weight  and  muzzle 
velocity,  the  latter  being  obtained  by  smear  camera.  Meteorological 
measurements,  comprising  winds  and  atmospheric  conditions,  were  taken 
throughout  the  time  of  the  firing. 


III.  REDUCTION  OF  FLIGHT  TEST  DATA 


At  Wallops  Island,  radar  data  consisting  of  slant  range,  elevation 
angle  and  azimuth  angle  are  recorded  on  digital  tape,  The  first  step  in  the 
actual  radar  reduction  program  is  the  application  of  the  appropriate 
trigonometric  transformations  to  convert  the  data  to  range,  height  and 
deflection  referenced  to  the  gun  coordinates  and  firing  azimuth.  These 
values  are  then  smoothed,  employing  tor  the  purpose  a  moving  arc 
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technique  that  uses  the  method  of  least  squares  to  fit  a  second-degree 
polynomial  to  the  data  in  every  7.  5  seconds  span  along  the  trajectory.  The 
polynomial  is  evaluated  at  the  midpoint  of  each  section,  and  the  first  and 
second  derivatives  are  computed  at  this  time  to  obtain  the  velocity  (u) 
and  acceleration  (u)  vectors. 


IV.  COMPUTATION  OF  DRAG  COEFFICIENT  (CD)  ' 


The  frame  of  reference  for  all  vectors  to  be  presented  is  a  ground- 
fixed,  right-handed  Cartesian  coordinate  system  with  unit  vectors 

•4  «4  <4  — $  *“♦ 

(1,2,3).  Define  the  2  axis  to  be  parallel  to  the  vector  g  at  the  origin 
and  in  the  opposite  direction. 

The  mass  of  the  projectile,  the  atmospheric  conditions  existing  at  the 

time  of  the  firing  a~  tn-  reduced  radar  data  for  the  round  provide  the 

necessary  inputs  to  the  equation  for  determination  of  the  drag  coefficient 

(C  ).  C  is  obtained  by  an  inverse  solution  of  the  point-mass  differential 
D  D 

equations  of  motion. 

CD  =  -l  (ttj  -  «,)(»!  -  -  At)  +  “2  (u2  -  g2  -  A2)  +  (u3  -  w3> 

2  3 

(u3  -  g  -  A3)J  8m/irpd  v 

-4 

The  acceleration  due  to  gravity  (g)  is 


*  -  <-fss>R2A3)  X2  +  R 


where  r  =  [  X  2  +  (Xj  +  R)2  +  Xj2]  *. 

The  Coriolis  acceleration  (A)  due  to  rotation  of  the  earth  is  expressed  by 

Fvz  •  v3] 


A  =  XjU,  +  X3u3 
^2U1  **  ^3U2 
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For  the  northern  hemisphere  the  A*s  are  defined  by  the  following  equations 
(for  the  southern  hemisphere  replace  L  by  -L). 


Aj  =  2  O  cos  L  sin  Az 
A,  -  2  f!  sin  L 

=  2  Cl  cos  L  cos  Az 


V.  RESULTS 


For  the  purpose  of  this  report  the  FPS  16  was  the  only  radar  that 
produced  acceptable  results  over  the  Mach  number  region  of  interest. 
Spandar  locked-on  very  late  in  the  flights  it  tracked,  yielding  too  few  data 
for  a  complete  analysis.  The  MPS  19  and  MOD  II  which  are  the  least 
precise  radars*'  gave  inadequate  results  using  this  method  of  analysis. 

All  of  the  data  collected  by  the  FPS  16  radar  were  reduced  and  used  in 
the  computation  of  aerodynamic  drag.  The  resultant  drag  coefficients  have 
been  grouped  according  to  various  combinations  of  quadrant  elevation  and 
charge  and  are  presented  in  graphical  form  as  a  function  of  Mach  number. 
Values  lying  outside  the  limits  of  any  one  graph  are  plotted  along  the 
edges  and  bear  the  following  identification  CD.  Figure  1  is  a  plot  of  the 
drag  coefficients  obtained  for  all  charges  and  quadrant  elevations.  As 
shown  by  the  graph,  the  curve  formed  by  these  computed  values  closely 
follows  the  aerodynamic  spark  range  curve  representing  the  expected 
performance  of  the  Ml 07  round. 

The  increase  in  drag  noticeable  in  the  lower  Mach  number  region  of 
the  plot  is  probably  due  either  to  a  combination  of  the  yaw  of  repose  and 
the  processional  and  nutational  yawing  motion  or  to  the  yaw  of  repose 
alone.  Figure  2  shows  the  subsonic  values  of  the  drag  coefficient  for 
Mach  numbers  less  than  .  8.  The  drag  coefficient  is  plotted  against 
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computed  values  of  the  yaw  of  repose  squared  as  defined  in  the  modified 

3 

point  mass  trajectory  model.  To  determine  the  change  in  drag  due  to 
yaw  a  least  squares  procedure  was  used  to  fit  these  data  as  a  linear 
function  of  yaw  squared.  The  resulting  curve  also  shows  good  agreement 
with  the  aerodynamic  spark  range  curve. 

Figures  3,  4  and  5  are  plots  of  the  aerodynamic  drag  coefficients 
found  at  each  of  the  quadrant  elevations  fired,  irrespective  of  charge. 
Figures  3  and  4  appear  to  indicate  that  at  the  lower  elevations  there  is 
little  if  any  increase  in  drag  due  to  the  yaw  of  repose,  while  at  the  highest 
elevation,  Figure  5,  the  increase  is  considerable.  The  remaining  Figures 
6  through  12  present  the  drag  coefficients  at  each  elevation  and  charge. 
Figures  6  and  7  re  plots  for  a  quadrant  elevation  of  300  mils  at  charges 
5  and  7.  In  Figures  8,  9  and  10,  quadrant  elevation  is  700 mils  at  charges  2, 
5  and  7.  Finally,  Figures  11  and  12  have  a  quadrant  elevation  of  1120  mils 
at  charges  5  and  7.  The  three  plots  at  charge  7  are  in  good  agreement,  and 
the  three  plots  at  charge  5  with  only  the  highest  quadrant  elevation  having 
a  large  amount  of  data  also  indicate  good  agreement. 

Figures  13  and  14  are  plots  of  the  drag  coefficients  calculated  for  all 
charges  and  quadrant  elevations  for  the  VPS  19  and  MOD  II  radars, 
respectively.  Both  graphs  indicate  that  meaningful  results  cannot  be 
obtained  from  these  radar  data  by  the  method  discussed  in  this  report. 


15 


VI.  CONCLUSION 


On  the  basis  of  the  above  study  it  has  been  established  that  it  i&  feasible 
by  a  polynomial  smoothing  and  differentiation  method  to  obtain  aerodynamic 
drag  force  values  from  flight  test  data  recorded  by  a  radar  with  the  pre¬ 
cision  of  the  FPS  16.  Assuming  a  computed  yaw  of  repose  as  defined  in  the 

3 

modified  point  mass  trajectory  model  ,  the  aerodynamic  change  in  drag 
due  to  yaw  can  also  be  determined.  Also,  it  can  be  concluded  that,  on  the 
average,  the  drag  results  obtained  from  radar  data  agree  well  with  the 
drag  results  from  ballistic  range  firings. 
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Yaw  Squared 

Figure  2,  Subsonic  Drag  Coefficient  (C^)  vs  Yaw  Squared 

for  All  Charges  and  QE's 
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Figure  4.  Aerodynamic  Drag  Coefficient  fC^)  vs  Mach  No 
for  Charges  2,  5  and  7,  QE  s  700  mils 


FPS  16 


FPS  16 


FPS  16 


REFERENCES 


1.  Wallops  Station,  Radar  Information  Handbook,  NASA. 

2.  R.  F.  Lieske  and  J.  W.  Kochenderfer,  "Determination  of  Performance 
Parameters  for  Fin-Stabilised  Free-Flight  Missiles,  "  Ballistic 
Research  Laboratories  Report  No.  1349,  AD  809790,  December  1966. 

3.  R.  F.  Lieske  and  M.  L.  Reiter,  "Equations  of  Motion  for  a  Modified 
Point  Mass  Trajectory,"  Ballistic  Research  Laboratories  Report 
No.  1314,  AD  485869,  March  1966. 


